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Abstract. We present 12 CO(1-0) and 12 CO(2-l) maps of the interacting Seyfert 2/LINER galaxy NGC5953 obtained with 
the IRAM interferometer at resolutions of 2'.'1 x 1'.'4 and l'.'l x C/.'T, respectively. We also present single-dish IRAM 30m 
observations of the central region of NGC 5953 for the 12 CO(1-0), 12 CO(2-l), and HCN(l-O) transitions at resolutions of 22", 
12", and 29", respectively. The CO emission is distributed over a disk of diameter ~16"(~2.2kpc), within which are several, 
randomly distributed peaks. The strongest peak does not coincide with the nucleus, but is instead offset from the center, ~2-3" 
(~340pc) toward the west/southwest. The kinematics of the molecular component are quite regular, as is typical of a rotating 
disk. We also compared the 12 CO distribution of NGC 5953 with observations at other wavelengths in order to study correlations 
between different tracers of the interstellar medium. The HST/F606W WFPC2 images show flocculent spiral structures and an 
"S-shape" feature i 60 pc in radius, possibly associated with a nuclear bar or with the radio jet. A two-dimensional bulge/disk 
decomposition of the /f-band (HST/F160W) and 3.6yum (Spitzer/IRAC) images reveals a circumnuclear "ring" ~10-14" in 
diameter, roughly coincident in size with the CO disk and with a star-forming ring previously identified in ionized gas. This 
ring is not present in the near-infrared (NIR) J - K color image, nor is it present in the "dust-only" image constructed from 
the 8 /um IRAC map. The implication is that the excess residual ring is stellar, with colors similar to the surrounding disk. 
We interpret this ring, visible in ionized gas, which appears as stars in the NIR, and with no sign of hot dust, as due to a 
red super giant population at least 10-15 Myr old. However, star formation is still ongoing in the disk and in the ring itself. 
Using NIR images, we computed the gravity torques exerted by the stellar potential on the gas. The torques are predominantly 
positive in both 12 CO(1-0) and 12 CO(2-l), suggesting that gas is not flowing into the center, and less than 5% of the gas angular 
momentum is exchanged in each rotation. This comes from the regular and almost axisymmetric total mass and gas distributions 
in the center of the galaxy. In NGC 5953, the AGN is apparently not being actively fueled in the current epoch. 

Key words, galaxies: individual: NGC 5953 - galaxies: spiral - galaxies: active - galaxies: nuclei - galaxies: ISM - galaxies: 
kinematics and dynamics 



1. Introduction 

Since molecular gas is the predominant phase of the interstellar 
medium (ISM) in the inner regions of spiral galaxies, CO lines 
represent an optimum tracer of nuclear gas dynamics, active 
galactic nuclei (AGN) fueling mechanisms, and their link with 
circumnuclear star formation. Although most galaxies host su- 
per massive black holes (SMBHs) and the gas accretion phe- 
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nomenon is usually invoked to explain nuclear activity in galax- 
ies, the nature of this activity is still not well known. The feed- 
ing of an AGN through accretion depends on an adequate sup- 
ply of gas whose angular momentum has been reduced enough 
to enable inflow on the small spatial scales surrounding the BH. 
Angular mo mentum must be removed from the disk gas (e.g.. 



Jogee 2006, and references therein), a process that can be ac- 



complished through non-axisymmetric perturbations of inter- 
nal or external origin. In the first case, they arise from disk 
instabilities and, in the se cond, from galaxy co llisions, merg- 
ers, and mass accretion dHeckman et al.l 1 1986b . Either way, 
they usually manifest themselves as density waves, such as 
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larg e-scale spirals or as bars and t heir gravity torques (e.g., 



Sakamoto e t al. 1999; Combes 2001), or as more localized phe 



nome na, including nested nuclear bars (e.g., Friedli & Martinet 



1993), lopsided ness or m — 1 perturbations (e.g., 



Shu et al 



1990; Garcfa-Burillo et al. 2000), or warped nuclear disks (e.g.. 
Schinnerer et al.ll2000allbl) 



To better understand the mechanisms for gas fueling of 
AGN, we started a high-resolution and high-sensitivity CO 
survey of nearby active galaxies at the IRAM Plateau de 
Bure In terferometer (PdBI) the N Uclei of GAlaxies (NUGA) 
project dGarcfa-Burillo et al.ll2003l) . The galaxies of the NUGA 
sample already studied show surprising results. In fact, there 
is no unique circumnuclear molecular gas feature linked with 
nuclear activity, but instead a variety of molecular gas mor- 
phologies which characterize the inner kpc of active galax- 
ies. T hese morphologies include one- and two-armed instabil- 
ities dGarcfa-Burillo et al.1 120031). well-ordered rings an d nu- 
clear spirals dCombes et al] l2004t Casasola et al] 2008), cir- 
cumn uclear asymmetrie s "dKrips et al.l 120051) and large-scale 
bars (lBoone et al.ll2007t Hunt et alj|2008l) . The analysis of the 
torques exerted by the stellar gravitational potential on the 
molecular gas of the NUGA sample has shown that the gas can 
be driven away from the AGN (e.g. for NGC 4321) or toward 
it (e.g. for NGC 2782, NGC 3147, and NGC 4579). However, 
the velocities observed for NUGA are too small to correspond 
to the AGN feedback models where violent molecular out- 



flows and superwinds are ex pected (e.g lNarayanan et aDl2006t 
Hopkins & Hernquistll2006l) . 



The different morphologies we find are probably re 



lated to the various timescales ( Garcfa-Burillo et aO 120051) . 



Large-scale bars can t ransport gas inward efficiently (e.g., 
ICombes & Gerinlll985t ISakamoto et al.ll 1999b . and it appears 



that they can also drive p owerful starbursts (e.g jKnapen et al 



2002; Jo gee et al. I l2005l ). Nevertheless, a clear correlation be- 



tween large -scale bars and nuclear ac tivity has not yet been 
found (e.g., Mulchaey & Regan] 1997 ). probably because the 
timescales for bar-induced gas inflow and AGN duty cy- 
cles are ver y different. Bars drive inflow over timescales 
(>300 Myr, Ijogee et al] 120051) larger than those of AGN 



accretion-rate duty cycles (~1-10 Myr, |H eckman et al. 2004 



iHopkins & HernquistH2006l iKing & PringleH2007l) . and active 
accretion seems t o occur only intermi t tently over the life 



2001: 



Mareckiet all [20031 



time of a galaxy (Ferrarese et al 

Janiuk et al ]|2004tlHopkins & Hernquistl2006t IKing & Pringl 
20071) . This implies that most AGN are in an intermediate 



phase between active accretion episodes, making the detection 
of galaxies with nuclear accretion rather difficult. 

Viscosity rather than self-gravity can also play a significant 
role in the fueling process. Viscous torques, generally weak 
and with timescales quite long at large radii, in combination 
with gravitational torques can coordinate efforts to produce re- 
current episodes of activity during the typical lifetime of any 
galaxy (IGarcfa-Burillo et al] 120051) . Viscous torques can pro- 
duce gas inflow on scales ~ 100-200 pc if they act on a con- 
trasted nuclear ring distribution and in the absence of partic- 
ularly strong positive gravitational torques. In NGC 4579, the 



efficiency of viscosity may be co mparable to the efficiency o f 
gravity torques in the inner ~50pc (IGarcfa-Burillo et alJ20 09). 

This paper, dedicated to the galaxy NGC 5953, is the latest 
of the NUGA series where results obtained for the galaxies of 
the sample are described on a case-by-case basis. NGC 5953 (D 
28 M pc for Hq _ 73 km s^Mpc -1 ) is an interac t ing galaxy 



(e.g., iRampazzo et al] 1 19951: [Casasola et all 12004 llono et al] 
2005), classified as a Seyfert 2 by Gonzalez Delgad o & Perezl 
d 19961) and as a LINER by IVeilleux et al.1 d 1995b . and of 
early and unbarred Hubble type (SAa pec). NGC 5953 and its 
late-type companion NGC 5954 (LINER/Seyfert 2, SAB(rs)cd 
pec) constitute a binary system (VV244, Arp91) where the 
two galaxies are separated by a projected distance of 5.8 
kpc (~43")- They show clear signs of interaction visible in 
the distorted morphology, the presence of a tidal bridge (or 
distorted arm) connecting the two galaxies, and of promi- 
nent star-forming regions. Both galaxies have circumnuclear 
starbursts that may have been induced by the interaction 
dGonzalez Delgado & Perezlll996l) . 

NGC 5953 hosts a compact radio core and jet, revealed by 
high-resolut i on rad io continuum observations with MERLIN 
( Kri ps et al. 2007al) . The jet is resolved at 18 cm, and af- 
ter beam deconvolution is roughly ~(J.'3 in length (40 pc), 
with a position angle (PA) of ~10°. The small-scale ra- 
dio continuum structure at 20 cm is similar in orientation, 
but with lower spatial resolution (~ l'/5), slightly more to- 
ward the east, PA~25° djenkinsl[l984]) . This jet-like elonga- 
tion as seen at lower resolution was fi rst referred to as a "jet" 
by iGonzalez Delgado & Perez (1996) because it is roughly 
aligned with the structure in the "excitation map", obtained by 
dividing [Om] emission by Ha. 

The NGC 5953/54 pair has also been mapped in atomic and 
molecular gas. Both galaxies are embedded in a common Hi 
envelope, with a clear velocity gradient alo ng the Hi plume 



extending more than 8 kpc to the northwest dChengalur et al 



1994: llono etaDl2005l;lHaan et al]l2007l l2008). There is some 



indicatio n of a faint diffuse optical counterpart to the H i 
plume dChengalur et al.lll994 [Hernandez-Toledo et alj|2003l) . 
The Hi peak is also s ignificantly displaced from the stellar 
disks dlono et al.ll2.QQ5h . The overall Hi velocity gradient runs 
from southeast to northwest, roughly perpen dicular to the rota 



tion in the ionized and molecular gas dHernandez-Toledo et al 



20031 llono et al.ll2005l) . The most recent Hi mass determina- 



tion for NGC 5953 has been obtained by lHaan et all d2008l) . 
Mhi = l.lxl0 9 M Q (value scaled to our adopted distance of 
D = 28 Mpc), typical of that exp ected for interacting galaxies 
of the same morphological type dCasasola et al.l l2004). 

The molecular gas distribution in NGC 5953 is symmet- 
ric with a shor t extension pointi ng toward NGC 5954 (e.g., 



Yao et alJ |2003; Iono et al ] [ 20051) . The H2 mass content es- 



timated by llono et al. (2005) is 2.2 x 1O 9 M (scaled to the 
distance of D = 28 Mpc for NGC 5953), higher than the 
molecular hydrogen mass found for the interacting compan- 
ion NGC 5954 (Mh, = 1.1 X 1O 9 M ). Higher-order transitions 
of the CO molecule have also been detected in NGC 5953, in- 
cluding the 12 CO(3-2) line by lYao et al] d2003l) with the James 
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Table 1. Fundamental parameters for NGC 5953. 



Parameter 


Value" 


References 


«J2000' 1 


15"34'"32.36 ! 


(1) 


<5j2000 J 


15°11'37.'70 


(1) 


Q'dyn 


15*34 ra 32.38 s 


(2) 




15°11'37.'59 


(2) 


^sys.hel 


1990 kms" 1 


(1) 


RC3 Type 


SAa pec 


(3) 


Nuclear Activity 


S2/LINER 


(4) & (5) 


Inclination 


42° 


(1) 


Position Angle 


45° ± 1° 


(1) 


Distance 


28Mpc (1" = 136 pc) 


(3) 


Lb 


4.9 x 10 9 L 


(6) 


M H i 


1.1 x 10 9 M Q 


(7) 


M H2 


2.2 x 10 9 M Q 


(8) 


M dusl (60 and lOOyUm) 


2.7 x 10 6 M 


(6) 


Lfir 


1.4 x 10 10 L Q 


(9) 


a ( a J20oo, £12000) is the phase tracking center of 


our 12 CO observ 



tions, (ttdyn. #dyn) is the dynamical cent er derived from rad io ob- 
servations for the core of NGC 5953 bv lKrips et al] J2007ah . 
Luminosity and mass values extracted from the literature have 
been scaled to the dista nce of D = 2 8 Mpc. 



(1) This paper; (2) iKrips et aT] 
Extragalactic Database (NED); (4) 



2007aj); (3) NASA/IPAC 



faxtrag atact ic Database (JNfcL)); (4) iGonzalez Delgado & PereJ 
1 1990): (5) [Veilleux et al] dl995h : (6) ICasasola et al] J2004) ; (7) 



lHaan et alj J2008h ; (8) llono et alj fc005h ; (9) IRAS Catalog 



Clerk Maxwell Telescope (FWHM~15"), suggesting a high 
excitation of the carbon monoxide, indicative of dense and 
hot gas. Table Q] summarizes the fundamental characteristics 
of NGC 5953. 

The structure of this paper is as follows. In Sect. [2j we 
describe our new observations of NGC 5953 and the litera- 
ture data with which we compare them. In Sects. [3] and [4] we 
present the observational results, both single dish and inter- 
ferometric, describing morphology, excitation conditions, and 
kinematics of the molecular gas in the inner kpc of NGC 5953. 
Comparisons between 12 CO observations and those obtained at 
other wavelengths are given in Sect. [5] In Sect. [6] we describe 
the computation of the gravity torques derived from the stellar 
potential in the inner region of NGC 5953. Sect. [7] summarizes 
our results. 

We will a ssume a distance to NGC 5953 of D = 28 Mpc 
(HyperLeda^, IPaturel et al. 2003 ) and a Hubble constant Hq = 
73 kms" 1 Mpc -1 . This distance implies that 1" corresponds to 
136pc. 



2. Observations 

2.1. Interferometric observations 

We observed NGC 5953 with the IRAM PdBI in the ABCD 
configuration of the array between September 2003 and 
February 2004 in the I2 CO(1-0) (115 GHz) and the 12 CO(2-l) 
(230 GHz) rotational transitions. The PdBI receiver character- 
istics, the observing procedu res, and the image reconstru ction 
are the same as described in iGarcfa-Burillo et alj (120031) . The 
quasars 3C454.3 and 3C273 were used for bandpass and flux 
calibrations respectively, and the quasar 1546+027 was used 
for phase and amplitude calibrations. 

Data cubes with 512 x 512 pixels (0''20 pixel for 
12 CO(1-0) and 0'.' 14 pixel" 1 for 12 CO(2-l)) were created over 
a velocity interval of -242 kms" 1 to +242 kms" 1 in bins of 
5 kms" 1 . The images wer e reconstructed using the standard 
IRAM/GILDAS0 software dGuilloteau & Lucasll2000r) and re- 
stored with gaussian beams of dimensions 2'.'0 x 1'.'4 (PA = 
25°) at 115 GHz and l'.'l x 077 (PA = 22°) at 230 GHz. We 
used natural weighting to generate the 12 CO(1-0) maps and 
uniform weighting to generate 12 CO(2-l) maps. Such a pro- 
cedure maximizes the flux recovered in I2 CO(1-0) and op- 
timizes the spatial resolution in I2 CO(2-l). In the cleaned 
maps, the rms uncertainty cr in 5 kms" 1 width velocity chan- 
nels is 2.7 mJy beam" 1 and 6.0 mJy beam" 1 for the 12 CO(1-0) 
and 12 CO(2-l) lines, respectively. At a level of ~3cr no 3 mm 
(1mm) continuum was detected toward NGC 5953, down to 
an rms noise level of 0.32 mJy beam" 1 (0.48 mJy beam" 1 ). The 
conversion factors between intensity and brightness tempera- 
ture are 32 K (Jy beam" 1 )" 1 at 1 15 GHz and 28 K (Jy beam" 1 )" 1 
at 230 GHz. All velocities are referred to the systemic velocity 
V syS; hei = 1990 kms" 1 (see later Sect. [47Tb and (Aa,A<5) off- 
sets are relative to the phase tracking center of the observations 
(a/2000, 6/2000) = (15 /, 34" , 32.36 s , 15°11'37770). All maps pre- 
sented in this paper are centered on this position. The maps are 
not corrected for primary beam attenuation. 

2.2. Single dish observations and short spacing 
correction 

We performed IRAM 30 m telescope observations in a 5 x 5 
raster pattern with 7" spacing from 16 to 19 July 2002. We used 
4 SIS receivers to observe simultaneously at the frequencies of 
the 12 CO(1-0) (115 GHz), the 12 CO(2-l) (230GHz), and the 
HCN(l-O) (89 GHz) lines. The half power beam widths are 
22" for 12 CO(1-0), 12" for 12 CO(2-l), and 29" for HCN(1- 
0). The CO(2-l) line has been observed in dual-polarization. 
Typical system temperatures were ~270-390K at 115 GHz, 
-300-750 K at 230 GHz, and -120K at 89 GHz. Throughout 
this paper the line intensity scale is expressed in units of 
r m b, the beam-averaged radiation temperature. T m \, is related 
to 7^, the equivalent antenna temperature reported above the 
atmosphere, by r\ — T*JT m \, where 77 is the telescope main- 
beam efficiency. At 115 GHz 77 = 0.79, at 230 GHz 77 = 0.54, 



1 http://leda.univ-lyonl.fr 



2 http://www.iram.fr/IRAMFR/GILDAS/ 
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Fig. 1. 12 CO(1-0) velocity channel maps observed with the IRAM PdBI+30m in the nucleus of NGC 5953, with a spatial reso- 
lution of 2'.'1 x 1'.'4 (HPBW). The maps are centered on the phase tracking center of our observations assumed to be coincident 
with the dynamical center of the galaxy (see Sect. 14. II ). Velocity channels range from AV = -130kms~' to +130 km s~' in steps 
of 5kms~' relative to V sys ,hei = 1990 kms~' (see Sect. 14. j} . The contours run from -0.20mJybeam _1 to 120mJy beairT 1 with 
spacings of 20 mJy beairT 1 . 



and at 89 GHz 77 = 0.82. All observations were performed in 
"wobbler-switching" mode, with a minimum phase time for 
spectral line observations of 2 s and a maximum beam throw 
of 240". The pointing accuracy was ~3" rms. The single dish 
maps are centered on the phase tracking center of the interfer- 
ometric observations (see TableHJ. 

Single dish 12 CO observations were used to compute short 
spacings and complete the interferometric measurements. We 
combined 30 m and PdBI data using the SHORT-SPACE task 
available in the GILDAS software. The best compromise be- 
tween good angular resolution and complete restoration of 
the missing extended flux was found by varying the relative 
weights of 30 m and PdBI observations. These were chosen 



in order to obtain the same mean weights in the single dish 
data as in the interferometric data within a ring in the uv plane 
ranging from 1.25 D/A to 2.5 D/A (D = 15 m). The combined 
PdBI+30m data have produced maps with angular resolutions 
of 2'.' 1 x 1'.'4 at PA = 25° for the 12 CO( 1-0) and V! Ixff.'l at PA 
= 22° for the 12 CO(2-l). In the combined maps, the rms uncer- 
tainty cr in 5 km s~' width velocity channels is 2.5 mJy beam -1 
and 5.5 mJy beairT 1 for the 12 CO(1-0) and 12 CO(2-l) lines, 
respectively. All figures presented in this paper are made with 
short-spacing-corrected data. 

We estimate that the 12 CO(1-0) map including only PdBI 
observations within -22" ( 12 CO(1-0) HPBW for the 30 m tele- 
scope) recovers a flux of Sco(i-O) = 150 Jy kms -1 , -81% of 
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Fig. 2. Same as Fig.[T]but for the 12 CO(2-l) line, with a spatial resolution of l'.'l x 077. The contours run from -60mJybeam 1 
to 150mJy beairT 1 with spacings of 30mJybeam _1 . 



the total flux measured with the combined PdBI+30m map, 
Sco(i-O) = 185 Jy kms- 1 . Table collects 12 CO(1-0) flux val- 
ues, both present in literature and determined with our obser- 
vations, for NGC 5953. In this table, Col. (1) indicates the ref- 
erence, Cols. (2) and (3) the telescope, both single dish and 
interferometer, and the diameter of the single dish telescope, 
Col. (4) is the primary beam of the instrument or the diam- 
eter used for the performed photometry, Col. (5) is the beam 
in interferometric measurements, and Col. (6) is the measured 
flux. Table [2] shows that 12 CO(1-0) fluxes we obtained with 
interferometric observations, single dish, and combined mea- 
surements (PdBI+30m) are in good mutual agreement with 
each other and with data present in literature. Our 12 CO(l- 
0) combined map within 42" gives a value (254 Jykms -1 ) 
consistent with t he flux of 233 Jy kms 1 found with OVRO 
dlono et alj|2005l) . but consistent also with the NRAO flux of 



365 Jykms" 1 dZhu et all 1 19991) taking into account that they 
used a single dish with a diameter of 12 m and a primary beam 
of 55". I n addition, we reco vered ~79% of the total flux de- 
tected by lYoung et aD d 19951) with the FCRAO (320 Jy km s 1 ), 
a reasonable agreement considering the uncertainties in the am- 
plitude calibration and the non-correction by the primary beam 
attenuation. 



Figures Q] and [2] show the channel maps for the 12 CO(l- 
0) and 12 CO(2-l) lines, respectively, in the central region of 
NGC 5953. Figures |4] and |5] display the single dish data, for 
the 12 CO(1-0), 12 CO(2-l), and HCN(l-O) lines. The two 12 CO 
lines have been mapped on a 5 x 5 grid with 7" spacings, while 
the HCN(l-O) line has been mapped on a 3 x 3 grid with 7" 
spacings and the nine HCN(l-O) spectra have been averaged to 
improve the signal-to-noise. 
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Table 2. 12 CO(1-0) flux values for NGC 5953. 



Reference 


Telescope 


Diameter 


Primary beam or FOV a 


Beam 


Flux 






[m] 


["] 


["X"] 


[Jy km s~'] 


Iono et al. (2005) 


OVRO 




60 


4.4 x 3.6 


233 


Zhuetal. (1999) 


NRAO 


12 


55 




365 


Young et al. (1995) 


FCRAO 


14 


45 




320 


This paper 


PdBI+30m 




42 


2.1 x 1.4 


254 


This paper 


PdBI+30m 
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2.1 x 1.4 


185 


This paper 


PdBI 
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2.0 x 1.4 


150 


This paper 


30 m 


30 


22 (central position) 




167 c 


This paper 


30 m 


30 


22 (inner 50" x 50") 




377 d 



a Primary beam is considered for single dish observations, while field-of-view (FOV) for interferometric or combined (interferometric+single 
dish) ones. 

b The photometry has been performed within 22", the I2 CO(1-0) primary beam for the 30m telescope. 
c The 12 CO(1-0) recovered flux for the central position (0", 0"). 

d The 12 CO(1-0) recovered flux for inner 50" x 50", 5x5 mapping with 7"spacing (see Sect. 12. 2\ . 



2.3. Optical and infrared images 

We first acquired a broad-band optical image from the HST 
archivd^ of NGC 5953 obtained with the F606W filter (mean 
wavelength o f 5940 A). This image was first published by 



A J - K i mage was derived from NIR data published 



Malkan et al 



(119981) . in a survey of 256 of the nearest 
(z<0.035) Seyfert 1, Seyfert 2, and starburst galaxies. The im- 
age covers the inner ~2Q" x 20" with a pixel size of 0'.'045. 

From the HST archive, we also acquired the HST NICMOS 
F160W (.ff-band, 1.6jum) image of NG C 5953. This near- 
infrare d ( NIR) image was prese nted by iRegan & Mulchaev 
d 1999b and lHunt & Malkanl(l2004l) . and covers the inner ~19"x 
19" of the galaxy with a pixel size of 0'.'075. 

We also acquired infrared (IR) images obtained with the 
IRAC camera on Spitzer, available thanks to the project 
'Starburst Activity in Nearby Galaxies' (Principal Investigator: 
G. Rieke). The IRAC images, from 3.6, to 8 fim, were reduced 
with MOPEX ( Ma kovoz & Marleau1l2005l) which accounts for 
distortion and rotates to a fiducial coordinate system. They 
cover a large sky area (~320"x320") including both NGC 5953 
and NGC 5954. We imposed a pixel size of 1"20 for the fi- 
nal images, roughl y the sam e as the original IRAC detector. 
Following iHelou et alj d2004i) . a "dust-only" (non-stellar) im- 
age was derived from the 8 |/m image, with subtraction of the 
stellar component computed by scaling the 3.6 yum and 4.5 fim 
images and subtracting them. This image will be referred to as 
a dust-only image, and should be dominated by emission from 
Polycyclic Aromatic Hydrocarbons (PAHs) and perhaps some 
hot dust continuum emission. 



3 Based on observations made with the NASA/ESA Hubble 
Space Telescope, and obtained from the Hubble Legacy Archive, 
which is a collaboration between the Space Telescope Science 
Institute (STScI/NASA), the Space Telescope European Coordinating 
Facility (ST-ECF/ESA) and the Canadian Astronomy Data Centre 
(CADC/NRC/CSA). 



by iHunt et all d 19991) acquired with ARNICA (Arcetri Near- 
Infrared Camera) mounted on the Nordic Optical Telescope 
(NO10). ARNICA was an imaging camera for the NIR bands 
between 1.0 and 2.5 jum based on the Rockewell HgTeCd- 
array detector NICMOS 3 (256x 256 pixels), and designed and 
built by Arcetri Observatory (Firenze, Italy) for the Infrared 
Telescope at Gornergrat in Switzerland. The ARNICA/NOT 
NIR images (J, H, K) cover ~120" x 120", including the pair 
of interacting galaxies, NGC 5953 and NGC 5954, with a pixel 
size of a/546. 

We compared the surface brightness profiles of the ground- 
and space-based images by extracting elliptically averaged pro- 
files, centered on the brightness peaks. The position angle and 
ellipticity were allowed to vary in the ellipse fitting. These ra- 
dial profiles are shown in Figure [3j the dashed horizontal lines 
in the lower panels correspond to the adopted position angle 
and inclination (see Sect. 14. II ). The profiles agree quite well; in 
particular the ground-based //-band and HSTF160W show the 
same trend down to the nuclear regions where the ground-based 
atmospheric beam smearing leads to lower surface brightness. 
The //-band/3.6 /im color is relatively constant throughout the 
entire radial range shown, with slightly redder colors from ~5- 



3. Single dish results 

The observations performed with the A and B receivers of the 
IRAM 30 m telescope in the two 12 CO lines covered the inner 
~50", corresponding to the central ~6.8 kpc (in diameter) of 
the galaxy (Fig. |4). The 25 observed positions show that the 



4 The NOT is operated on the island of La Palma jointly by 
Denmark, Finland, Norway, and Sweden, in the Spanish Observatorio 
del Roque de los Muchachos of the Instituto de Astrofisica de 
Canarias. 
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Fig. 3. Radial surface brightness profiles of NGC 5953 made by elliptical averages. The dotted line corresponds to IRAC 3.6 /im, 
the dashed line to ground-based H-band, and the solid line to HST/NICMOS/F160W. The left panel shows a smaller FOV than 
the right panel. The adopted position angle (45°) and the best-fit inclination (42°) are shown by dashed horizontal lines in the 
lower panels (see Sect. 14. U . The dashed vertical lines correspond to the 60 pc (~(y.'44, radius) nuclear bar seen in the HSTF606W 
image (see Sect. 15. U . and the 820 pc (radius) ring observed in the residuals of the bulge-disk decomposition (see Sect. 15.2b . 



central region of NGC 5953 hosts extended molecular emis- 
sion both in 12 CO(1-0) and 12 CO(2-l) emission (Fig. Hi. The 
maximum detected T m t is 0.18 K in 12 CO(1-0) in the south- 
west corresponding to offset position (-7", -7"), and 0.54 K in 
12 CO(2-l) at the same offset position (-7", -7"). 

Within the inner ~50" x 50" we estimate a flux of 377 
Jy km s _1 , in g ood agreement with previous singl e dish flux de- 
terminations dZhu et al. I ll 99911 Young et all 19951 see Table H. 
From this 12 CO(1-0) integrated flux assuming a H2-CO con- 
yersion factor X = N(U 7 )/ I C n = 2.2 x 10 20 cnr 2 (K km s" 1 )" 1 
dSolomon & Barrettlll99lh . we can derive the H2 mass within 
the observed region as: 



M H JM ] = 8.653 x 10 3 D 2 [Mpc]5 CO (i-o)[Jykms- 1 ] 



(1) 



We obtain Mh,~2.6x1O 9 M , and including the mass of helium, 
the corresponding total molecular mass is M mo i = Mn 2+ He = 
1.36xM H2 ~3.5xl0 9 M G . 

The HCN(l-O) line has been observed for 9 positions with 
7" spacing, covering the central ~43" (~5.8 kpc). The HCN(1- 
0) emission is detected in the west part of the observed region, 
and the average spectrum over the 3x3 grid is shown in Fig. [5] 
The HCN(l-O) peak is at r mb ~0.005 K. The CO(1-0)/HCN(1- 
0) ratio is very high, equal to 20 on average over the center, 
where instead we would expect an enhanced HCN emission 
due to the AGN and thus a lower ratio. A high CO/HCN line 
ratio is rarely encountered in AGN but not completely unusual. 



NGC 695 1 has a ratio of 30 in the starburst ring and 2.5 in the 
nucleus dKrips et al.ll2007bh . and NGC 3147 of 20 in the inner 
4 kpc where the CO emission exhibits tw o ring-like structures 
around the nucleus (ICasasola et al.ll2.Q08b . On the other hand, 
the galaxy NGC 1097 has a CO/HCN ratio which ranges from 



3 in the nucleus to 10 in the star-forming ring (Kohnoetal 



20031) . Since AGN activity implies an enhanced HCN emission 
relative to CO emission associated with the star-formation pro- 
cess, the high CO(1-0)/HCN(1-0) ratio observed in NGC 5953 
implies that excitation by star-formation is dominant over AGN 
excitation in the circumnuclear region. 

4. Interferometric results 

4.1. Dynamical center and inclination 

We assume that the dynamical center of NGC 5953 coin- 
cides with the position of the AGN "c ore" derived from ra - 
dio observations (MERLIN/ 18 cm) by iKrips etall d2007ah : 
ff=15 /, 34'"32.38 s and <5=15°1 1'37'.'59. Since these coordinates 
are nearly coincident with those of the phase tracking center of 
our 12 CO observations (see Table [TJ, in the following we as- 
sume that our observations are centered on the dynamical cen- 
ter of the galaxy. 

The spectral correlators were centered at 114.520 GHz 
and 229.037 GHz for the 12 CO(1-0) and 12 CO(2-l) line re- 
spectively, corresponding to Vlsr = 1966 kms~'. The differ- 
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Fig. 4. Spectra maps of NGC 5953 made with the IRAM 30 m with 7" spacing in 12 CO(1-0) (top) and 12 CO(2-l) (bottom). The 
positions are arcsec offsets relative to the phase tracking center of our observations. Each spectrum has a velocity scale from 
-300 to 300kms _I , and abeam-averaged radiation temperature scale (r m b) from -0.10 to 0.37 K for 12 CO(1-0) and from -0.25 
to 0.75 K for 12 CO(2-l). 



ence between LSR and heliocentric velocity is 16 kms , and 
therefore the observations were centered on Vhei(PdBI) = 1950 
km s . In the inner ~4" the velocity centroid is 50 km s red- 
shifted with respect to the heliocentric velocity of the center of 
our 12 CO(1-0) observations (Fig. [6] left panel, one component 
gaussian fit), and 30 kms -1 redshifted relative to 12 CO(2-l) 
(Fig- HI right panel, one component gaussian fit). We therefore 
estimate the systemic heliocentric velocity as the intermediate 
value between that determined for the 12 CO(1-0) and that for 
the 12 CO(2-l), V sys ,hei = 1990 kms" 1 . 



In Figs. [6] we fitted the emissions also considering 
three gaussian components, visible both in 12 CO(1-0) and in 
12 CO(2-l). These (red) fits show that the wing components 
for the two 12 CO lines are approximately at the same veloc- 
ities within the noise (-52 km s -1 for the 12 CO(1-0) vs. -57 
km s- 1 for the I2 CO(2-l), and 118 km s" 1 for the 12 CO(l- 
0) vs. 113 km s _1 for the I2 CO(2-l)). The major difference is 
present for the central gaussian component, 33 km s _1 for the 
12 CO(1-0) vs. 18 km s _! for the I2 CO(2-l), maybe due to the 
clumpy nature of the inner molecular gas. The line ratio R21 
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Fig. 5. HCN(l-O) spectrum toward the center of NGC 5953, 
averaged over the 9-point map made with the IRAM 30 m with 
7" spacing. The spectrum has a velocity scale from -500 to 
500 km s' 1 and a beam-averaged radiation temperature scale 
(r mb ) from -0.0035 to 0.006 K. 



(Jco(2-\)l lco(\-0)) assumes the values of 1.2, 0.6, and 0.5 re- 
spectively for the three components passing from negative to 
positive velocities. The mean ratio of 0.8 is consistent with the 
line ratio discussed later in Sect. [ 



The systemic heliocentric velocity of 1990 kms -1 is 25 
kms -1 redshifted with respect to the value determined from 
Hi observations (1965 kms 1 , HyperLeda), and in agreement 
with the shift of 50 kms - 1 inferred from p revious 12 CO(1-0) 
observations with OVRO dlono et al.ll2.005D . This discrepancy 
between systemic velocity derived from 12 CO and Hi obser- 
vations is not unusual, especially in galaxies with a lopsided 
Hi morphology. The NUGA galaxy NGC 4579 exhibits a dif- 
ference of ~50 kms -1 , probably due to the interaction his- 
tory of the galaxy and the differe nt effect of the ram-pressur e 
on the atomic and molecular gas ( Gar cfa-Burillo et al. 2009). 
Since the atomic gas is much more a ffected by ram-pressure 



(Kennev & Young 1986; Vollmeret al. 2001) than the molec 



ular component, the Hi kinematics are more sensitive than the 
CO to the interaction history; NGC 5953 could present a simi- 
lar scenario to NGC 4579. In fact, in NGC 5953 the Hi has been 
detected mainly in the outer parts where a big tidal tail connect- 
ing the two galaxies shifts the Hi barycenter with respect to the 
molecular gas that traces the less perturbed, inner parts of the 
galaxy. 

We also derived the inclination of NGC 5953 by fitting 
the two-dimensional bulge/disk (B/D) decomposition on the 
ground-based //-band image and the IRAC 3.6 /jm image (see 
Sect. 15.2b . The best-fit inclinations are i = 39° and 46°, respec- 
tively; we therefore used a roughly intermediate value of 42°, 
which approximates quite well the outer regions of the NIR 
images (see Fig. [3). 



4.2. CO morphology and mass 

The I2 CO(1-0) and I2 CO(2-l) integrated intensity distribu- 
tions are shown in Figure [7] The CO emission is distributed 
over a disk of ~16"(~2.2kpc) diameter. Our 12 CO(1-0) ob- 
servations show several peaks, distributed more or less ran- 
domly, with the strongest one offset from the nucleus ~2" to- 
ward the west/southwest, away from NGC 5954. The 12 CO(l- 

0) distribution is different from the Hi morphology, where the 
strongest emission has b een detec ted mostly on the side nearer 
to NGC 5954 (e.g Jlono et al.l2005l) . although the Hi resolution 
is probably insufficient to resolve distinct peaks. In the 12 CO(2- 

1) map the central emission is also clearly resolved and more 
clumpy than in I2 CO(1-0). The strongest I2 CO(2-l) peak is 
not that at ~2-3" in the west/southwest direction from the nu- 
cleus, like for I2 CO(1-0), but that at ~1"5 in the east direction 
from the nucleus. 

The total H2 mass derived from our interferometric 
(PdBI+30m) map within the 42" primary beam field of the 
PdBI (S co = 254 Jy kms -1 , see Table H is M H2 ~1.7xl0 9 M 
(M mo i~2.3xl0 9 Mq). Within a diameter of 22" (S C o = 185 Jy 
kms 1 , see Tabled, we derive a H2 mass of Mh,~1.3x10 9 M , 
~50% of the single dish 30 m H2 mass (see Sect. [3) corre- 
sponding to a region about 6.6 times larger in area. Hence, half 
of the molecular gas is concentrated in the central region of 
NGC 5953, as expected for spiral galaxies. The molecular gas 
distribution consists mainly of individual giant molecular cloud 
complexes, with the biggest one alone having a mass of Mh 2 = 
3.0 x 10 8 M Q , while the smaller ones have masses of a few 
10 7 M . 

NGC 5953 is quite massive in molecular gas compared to 
the other NUGA galaxies, where typically the mass is on the 
order of ~ 3xl0 8 Mq. Three g alaxies of th e NUGA sample , 
NGC 4569 (iBoone et al.ll2007h. NGC 2782 jHunt et allbooi) . 
NGC 3147 dCasasola etal.120 08). with molecular masses in the 
range of l-3xl0 9 M Q , are com parable with NGC 59 53. The ex- 
traordinary case of NGC 1961 dCombes et al.l2009l) is excluded 
from these considerations: NGC 1961 exhibits a H2 mass of 
~1.8xl0 10 Mq, almost an order of magnitude more massive 
than any NUGA galaxy. 



4.3. CO line ratio 

The comparison between the two 12 CO maps, obtained after 
convolving the I2 CO(2-l) map to the lower resolution of the 
12 CO(1-0) one and including short spacings for both maps, 
gives information about the local excitation conditions of the 
molecular gas. Figure [8] shows the ratio of the 12 CO(1-0) to 
I2 CO(2-l) convolved to the same resolution with 12 CO(1-0) 
contours as in Fig.|7](left panel). The line ratio ranges from 0.3 
to 0.9 inside the observed region and the bulk of the emission 
has a ratio between 0.4 and 0.8, value s consistent with optically 
1 s piral d isks (e.g.. iBraine & Combeslll992t 
ail l 19931) 



thick emission in 



Braine & Combeslll992 
Ratio values of ~0.9 are reached 



Garcfa-Burillo et 
in some clumpy structures, randomly distributed inside the 
mapped region, especially those toward E/NE with respect to 
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Fig. 7. Left panel: 12 CO(1-0) integrated intensity contours observed with the IRAM PdBI+30m toward the center of NGC 5953. 
The white star marks the coordinates of the dynamical center of the galaxy, with offsets in arcseconds. The map, derived with 
2<t clipping, has not been corrected for primary beam attenuation. The rms noise level is cr — 0.09 Jy beam -1 km s _1 and contour 
levels run from 3cr to 2\cr with 3cr spacing. In this map the ±130 km s velocity range is used. The beam of 2'.'1 x 1'.'4 is plotted 
in the lower left. Right panel: Same for 12 CO(2-l). The rms noise level is cr — 0.2 Jy beam 1 km s and contour levels run from 
3cr to lOcr with 3<r spacing. The beam of l'.'l x 077 is plotted at lower left. 



the nucleus, well visible in 12 CO(2-l) emission (Fig. [7] right 
panel). A ^21 line ratio of 0.9 suggests a mode rately high 
12 CO excitation consistently with lYao et al.l (120031) who have 
detected the 12 CO(3-2) line (345.80 GHz) with an intensity of 
hi = 38.7 K km s _1 and a line ratio R32 = /32//10-0. 85 in 
the nucleus of NGC 5953. Different excitation conditions of the 
molecular gas appear to characterize the interacting companion 



galaxy NGC 5954, where the 12 CO(3-2) line has been detected 
with a ratio of # 32 =0.37 dYao et alJl2003h . 

4.4. Kinematics 

The velocity-channel maps (see Figs. Q]and[2]i show a general 
regularity of the large scale kinematics, typical for a rotating 
disk. Some local wiggles, not forming a coherent grand design, 
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Fig. 8. Color scale of the CO(2-1)/CO(1-0) ratio map and 
12 CO(1-0) intensity map contours as in Fig.[7](left panel). 



are superimposed on this regular pattern, both at negative and 
positive velocities especially to the west of the nucleus. These 
kinematic "glitches" are probably associated with the intensity 
peak toward the SW. 

Figure [9] shows the 12 CO(1-0) isovelocity contours (first- 
moment map) superposed on the 12 CO(1-0) integrated inten- 
sity. The white star indicates the dynamical center of the galaxy 
assumed coincident with the phase tracking center of our ob- 
servations, and the velocities are relative to the systemic helio- 
centric velocity, V sys ,h e i = 1990 kms -1 . The dotted line traces 
the major axis of the galaxy determined from our observations 
(PA = 45° + 1°), by maximizing the symmetry in the position 
velocity diagrams. This value is also the best-fit PA given by 
the elliptically-averaged surface brightness profiles discussed 
in Sect. 1231 



Position-velocity (p-v) cuts along the major (PA = 45°) and 
minor axis (PA = 135°) of NGC 5953 are shown in Figures [10] 
andfTTl respectively. In both figures, the 12 CO(1-0) emission 
is given in the top panel and the 12 CO(2-l) one in the bot- 
tom panel. The p-v diagrams along the major axis reveal regu- 
lar circular rotation and appear symmetric in the inner 10-11", 
both in 12 CO(1-0) and 12 CO(2-l). This regularity in the p-v 
diagram of NGC 5953, more typical for normal spiral galaxies 
than interacting ones, has also been found in the ionized gas, 
includi ng Ha, [OIII1 JGonzalez Delgado & Perez! [r 996). and 



[Nil] dHernandez-Toledo et al J 12003). The kinematics shown 
by the minor-axis p-v diagrams (Fig. [TTJ are also quite reg- 
ular, but show significant velocity dispersion for the 12 CO(l- 
0) close to the nucleus. This effect can be attributed to beam 
smearing. 



Fig. 9. Overlay of the integrated 12 CO(1-0) emission, same as 
Fig- 13 (left panel), with CO mean-velocity field in contours 
spanning the range -100 to 100 km s _1 in steps of 10 km s -1 . 
The white star indicates the dynamical center of the galaxy. The 
velocities are referred to V sys .hei = 1990 km s _1 , solid (red) lines 
are used for positive velocities, and dashed (blue) lines for neg- 
ative velocities. The dotted line indicates the position angle of 
the major axis (PA = 45° + 1°). 



4.5. The rotation curve and dynamical mass 

We have derived a rotation curve (RC) from the p-v diagram 
taken along the kinematic major axis of NGC 5953 at a PA — 
45 °. The terminal velocities were derived by fitting multiple 
gaussian profiles to the spectra across the major axis. The fitted 
velocity centroids, corrected for sin i (i = 42°), give V Q \, S / sin i 
for each galactocentric distance. 

For both lines the velocities (negative) to the northeast in- 
crease less steeply within ~250pc than those to the south (pos- 
itive velocities). Nevertheless, they converge to the same value 
at greater distances, so we averaged together the two curves de- 
rived from either side of the major axis. Because data for both 
lines were consistent, we combined both curves into an average 
by spline interpolation. The RCs from the two 12 CO transitions 
are shown in Figure [T2] They are consistent with those found 
for the ion ized gas, as reported bvlponzalez Delgado & Perez 
1996) and Hernan dez^Toledo et all d2003l) . 



The 12 CO RCs are very regular and apparently this be- 
havior is expected more for normal spiral galaxies than inter- 
acting systems. NGC 5953 and the companion are clearly in- 
teracting and their outer parts are perturbed as shown by the 
presence of a tidal tail seen both in op t ical and in Hi em is- 
sion (e.g. jHernandez-Toledo et af 2003 : Haan et al. 2008 ). In 
NGC 5953, we are mapping the molecular gas in the inner parts 
that have a dynamical time-scale much shorter than the outer 
ones, and therefore they had the time to relax and reach an equi- 
librium state, perhaps after gas accretion/exchanges from/with 
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NGC 5954. In other words, the regularity of the RCs in 12 CO is 
not unexpected. 
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Fig. 10. Upper panel: I2 CO(1-0) position-velocity diagram 
along the major axis (PA = 45°) of NGC 5953 using the veloc- 
ity range from -140 to 140 km s contoured over a grey-scale 
representation. Contour levels are from 3<x to 39cr in steps of 
2cr (cr=2.5mJybeairr'). The velocities are relative to V S y S ,hei 
(=1990 kms -1 ) and X are the offsets along the major axis in 
arcsecs. Bottom panel: The same for 12 CO(2-l). Contour lev- 
els are from 3<x to 20cr in steps of 2<x (cr=5.5 mJy beam -1 ). 



In the major-axis p-v diagrams, the velocity peaks at 
~120/(sin/)kms _1 at a radius of ~2-A" (~400pc); however, 
the bulk of the gas is rotating at ~98 / (sin z) km s (see Fig.fT2li. 
This is roughly consistent with what is found for the ionize d 

" |2003fl 



gas, ~115/(sin/)kms by Herna ndez-Toledo et al 



The bulk velocity corrected for disk inclination, sin i (we adopt 
i = 42°, see Sect. I4.lt . becomes ~146kms~ 1 . From the cor- 
rected bulk velocity we can estimate the dynamical mass within 
a certain radius using the formula M(R) = axRV 2 (R)/G, where 
M(R) is in M , R in kpc, and V in kms~'. Assuming the most 
flattened disk-like distribution (a = 0.6), the above formula 
gives a dynamical mass of Md y „ = 2.4 x 10 9 M Q within a radius 
of 6" (-0.8 kpc). 
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Fig. 11. Upper panel: Same as Fig[10]along the minor axis (PA 
= 135°) of NGC 5953. Contour levels are from 3cr to 22cr in 
steps of 2<x. Y are the offsets along the minor axis in arc- 
secs. Bottom panel: The same for 12 CO(2-l). Contour levels 
are from 3<x to 12<x in steps of 2<r. 



Assuming a roughly flat rotation curve, the dynamical 
mass should be Md yn = 8.6 x 1O 9 M within a radius of 21" 
(~2.9kpc). In the same region, we estimate a total molecular 
gas mass of ~1.7xl0 9 M o (see Sect. l4~2li. a fraction of -20% 
of the dynamical mass. Anyway, this percentage is subject to 
uncertainties. One of these could be due to the variation of the 
H2-CO conversion factor. In this work, as for similar studies 
of molecular gas, we assume a constant value for the X con- 
version factor, for instance that determined for the Milky Way. 
However, the X ratio is determined by various factors, such as 
the metallicity, the temp erature, the cosmic ray density, and the 
UV r adiation field (see iMaloney & Blacklll988l: iBoselli et alJ 



2002), and its value can cha nge by a factor 4-15 (e.g., Wilson 



Taking into account their different adopted inclination. 



1995c ICasasola et alJ 120071) . The X value also changes with 
galaxy morphological type: usually galaxies earlier than Scd 
show values comparable to, or lower than, the Galactic one, 
while extremely late- type spirals or irregu lar galaxies tend to 
show higher values dNakai & Kunol 1995 ). Another possible 
uncertainty, but on the estimate of the dynamical mass, could 
be the suspected warp of the disk, discussed later in Sect. 15.31 
that might give a wrong inclination of the disk, and so an in- 
correct rotation velocity. 
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Fig. 12. Left panel: The CO(l-O) rotation curve derived as described in the text. The positive (negative) velocities are shown as 
a red dotted (blue dashed) line; the open (green) circles show the average. The horizontal dotted line at 146 km s indicates the 
velocity of the bulk of the molecular gas. Right panel: The same as the left panel, but for 12 CO(2-l). 



5. Dust, stellar structure, and star formation 

In this section we compare the 12 CO distribution of NGC 5953 
with observations at other wavelengths. These comparisons 
have a dual aim. First, they allow the study of possible cor- 
relations between different tracers of the ISM, such as between 
molecular clouds, considered the birth site of future genera- 
tions of stars, and the light from the stellar component. The sec- 
ond aim is to determine the location of dynamical resonances, 
which greatly aid the determination of gas flow in the circum- 
nuclear regions of the disk. However, in the case of NGC 5953, 
we find no evidence for any non-asymmetric component which 
could drive resonances. 



5.1. Optical and near-infrared morphology 

Figure [O] (left panel) displays a HST/NICMOS/F606W band 
image of the inner 20" of NGC 5953, described in Sect. 12.31 
We imposed that the brightness peak in this image coincides 
with the phase tracking center of our 12 CO observations and 
therefore with the dynamical center, as discussed in Sect. 14.11 
The small peak at ~3" west of center is al most certainly 



the for eground star identified in NGC 5953 by iRafanelli et al 



(1990). nspection of this figure shows a clearly-defined non- 
axisymmetric structure of ~250pc in size, similar in shape and 
size to the [Om] emission found by lGonzalez Delgado & Perez 
d 1996b . A close-up of the inner 2" (Fig.[T3] right panel) reveals 
an "S -shape" feature, perhaps assoc iated with a nuclear bar 
(e.g., iHernandez-Toledo et alJ 120031) . This bar-like signature 



is it 6 pc in radius, with a PA~17°. IHernandez-Toledo et al 



(2003), by examining a contrast-enhanced B-band image of 
NGC 5953 (the same image as used here), identified this struc- 
ture at about the same PA as the global major axis of the 
galaxy. We find (see Fig. [13] right panel) instead a more ver- 
tical orientation, PA~5 - 20°, depending on which features of 
the S-shape are used to define the PA. Rather than a nuclear 



bar, the S-shape stru cture could also b e an optical counter- 
part of the radio jet dKrips et alJl2007al) . because of the sim- 
ilar morphology, position angle, and size; this feature also 
seems to correspond to a higher e xcitation of the ionized gas 
dGonzalez Delgado & Perezll 19961) . 

The inner 1.6 pm morphology (HST/NICMOS/F160W) of 
NGC 5953 is shown in Figure[14] The foreground star to the 3" 
to the west of the nucleus is clearly visible also at 1 .6 pm. The 
fiocculent spiral structure clearly visible in the optical (Fig. [13} 
is not seen in the NIR emission; there is also no NIR counter- 
part to the S-shape feature seen in the optical. 

The morphology of the warm dust in the circumnuclear re- 
gions of NGC 5953 is shown in Figure [BJ where 12 CO(1-0) 
and 12 CO(2-l) intensities are contoured over the dust-only im- 
age described in Sect. 12. 31 There is a dust emission peak to the 
SW, not exactly coincident with the 12 CO emission peak but 
roughly in the same direction. The bulk of the (putative) dust 
emission is configured in a smooth featureless disk, similar to 
that seen in the 12 CO emission. 

5.2. Stellar structure 



A large-scale view of the interacting companions NGC 5953 
and NGC 5954 at 3.6pm is illustrated in Figure [TBI Like the 
K band, this image traces the massive component of the stellar 
populations in galaxies, but at 3.6 um the extinction is about 
1/3 that of the K band dCardelli et al.lll989b . 

To investigate stellar structure, we p erformed a two- 
dimensional B/D decomposition with galfit dPeng et alj|2002l) 
on both the ground-based //-band image and IRAC 3.6 pm im- 
age, fitting both galaxies (NGC 5953 and NGC 5954) simulta- 
neously. In the fit, the models were convolved with the IRAC 
Point Response Function (PRF) or an approximation of the 
ground-based Point Spread Function, as appropriate. The best 
fit included a Sersic generalized exponential bulge, an expo- 
nential disk, and a nuclear point source for NGC 5953, and a 
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Fig. 14. Left panel: F160W band NICMOS/HST image of NGC 5953. The inner 12" are shown. Right panel: A zoom of the 
nuclear region. The inner 3" are shown. 



Sersic generalized exponential bulge (but with index n ~1, so 
this is probably in reality a disk) and a nuclear point source 
for NGC 5954. The IRAC fit has a reduced^ of 4.2 (averaged 
over NGC 5953 and NGC 5954), roughly ten times smaller than 
that in the //-band. Hence, in what follows, we mainly rely on 
the IRAC fit. The bulge comprises more than half (~59%) of 
the total 3.6 pm light, consistently with its early Hubble type, 
but has a small (n ~ 1) Sersic index, unlike most early-type 
bulges with n ~ 3 dMoriondo et al . 1998). The disk scalelength 



is ~ 4kpc, slightly lar ger than typical early-type spiral disks 
(e.g.. lHuntetal]|2004 : this could be due to the interaction and 
the proximity with NGC 5954. 

The IRAC 3.6 pm residuals from the B/D decomposition 
are shown in Figure[17] with the 12 CO(1-0) (left panel) and the 
12 CO(2-l) (right panel) in contours. The residuals show a low- 
amplitude circumnuclear "ring" with a radius of ~5-7". The 
12 CO(1-0) and 12 CO(2-l) disks nearly coincide in size with 
the ring structure in the residuals. 
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Fig. 17. Left panel: I2 CO(1-0) integrated intensity contours as in Fig. [7] (left panel) overlaid on the residuals (in MJy sr _1 ) of 
the B/D decomposition performed on the IRAC 3.6 pm image. The (red) cross marks the dynamical center. The central 30" are 
shown. Right panel: Same for 12 CO(2-l) integrated intensity contours as in Fig. [7] (right panel). The central 20" are shown. 



Because the size of the ring in the IRAC residuals is typi- 
cal of small-scale "doughnuts" associated with the diffraction 
limit and PRF incompatibilities, we wanted to verify the re- 
ality of the ring-like feature seen in the 3.6 pm residuals. If 
real, it should also appear in the //-band residuals, which are 
shown in the left panel of Figure [18] Inspection of the figure 
shows that the ring is indeed present, delineated by a series of 



surface brightness excesses. The intensity contrast in these is 
around a factor of 10 relative to the regions external to the ring. 
The presence of the ring in both sets of residuals suggests that 
it is a real feature, rather than an artefact of IRAC PRF mis- 
matching. Moreover, intensity peaks in the ring-like structure 
are also present in the un-sharp masked HST/NICMOS/F160W 
map (see Hunt & Malkan 2004). We therefore conclude that the 
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Fig. 18. Left panel: //-band residuals from the galfit B/D decomposition. Right panel: The J - K color image. The FOV for both 
panels is 20" with north up, and east to the left. In both images, the foreground star ~2'.'7 to the west is evident. 
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Fig. 16. IRAC 3.6 pm image of NGC 5953 and NGC 5954, co- 
ordinate system centered on NGC 5953 (with N up, E left). The 
inner 120" are shown. 



ring in the NIR residuals is real, and probably stellar in nature, 
as we discuss below. 



5.3. The stellar ring 

The right panel of Fig. [18] displays the ARNICA/NOT J - K 
color image of NGC 5953. The nucleus of NGC 5953 and a 
small region to the west have extremely blue J - K colors 



(~0.7 and ~0.3, respectively), while the disk has rather nor- 
mal ones (~1.0). The small knot 3" west of center, visible also 
in the WFPC2/F606W (Fig. [T3 and NICMOS/F160W (Fig. 
[T4l i images, coinc i des al most certainly with the star identified 
bv lRafanelli et al. (1990). There may be a few isolated regions 
(~1") to the east of the nucleus with a red excess ( J - K~1.3), 
but there is no evidence in the J - K colors for the ring seen in 
the B/D residuals. This suggests that the excess residual ring is 
a stellar ring, with colors similar to the surrounding disk. 

A hint of the circumnuclear ring in NGC 5953 also appears 
in previous optical observations of the ionized gas. With the 
integral-field spectr ograph SAURON operating o n the William 
Herschel Telescope, iFalcon-Barroso et al] d2006 f ) found recent 
star formation in a circumnuclear ring with a radius of ~6- 
7". The H/3 emission is apparently distributed in a circum- 
nuclear disk, roughly coincident with the CO one, but the 
H/3 velocity dispersion and the [Ora]/HjS ratio are lower in a 
ring-like structure than in the main circumnuclear disk. This 
could be because the nucleu s excites gas in an extended region 
( Falcon-Barroso et alj|2006l) . or be due to a true ring in the dy- 
namical sense of a resonance. Lower velocity dispersion and 
lower [Oml/ Hff are both associated w i th enh anced star forma- 
tion, leading iHernandez-Toledo et al.1 (12003b to conclude that 
NGC 5953 hosts a circumnuclear ring. The size of this star- 
forming ring is coincident with that of the NIR ring in the resid- 
uals. Figure [19] shows the superposition of the 3.6/mi residu- 
als image (Fig. [17] left panel), represented in (green) contours, 
on the //-band residuals image (Fig. [18] left panel). This fig- 
ure clearly shows the coincidence in size between the two NIR 
rings, the same size of the star-forming ring seen in the ionized 
gas. 

Because the ring is visible in ionized gas, but in the 
NIR appears as stars, with no sign of hot dust because 
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presence of the Hi tidal tail gives some hint that the interaction 
in the pair NGC 5953/5954 is prograde, i.e., that at least one 
of the disks involved has its spin axis parallel to the orbital ve 



locity of the merger dHibbard & van Gorko m 1996; Ionoetal 
l2005h . 



5 0-5 
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Fig. 19. 3.6 //m residuals in (green) contours (Fig.[T7li overlaid 
on the //-band residuals (Fig. [18] left panel). The inner 20" are 
shown. 



of normal J - K colors, we interpret it as an ongoing 
starburst, with a stellar Red Super Giant (RSG) popula- 
tion. This would imply that the ring feature is at least 
10-15Myr old , beca use RSGs onset at roughly that time 
(Leit herer et al. 1999h . Ho wever, the presence of hydrogen 
recombination lines (e.g., Gonzalez Delgado & Perez 19961; 



Hernandez-Toledo et al 2003) means that star formation is still 
ongoing, although the exact age is difficult to determine (5- 
lOMyr). 

In a ddition to the star-forming+ stellar ring, the observa- 
tions of iFalcon-Barroso et al.l (12006b suggest that the stars in 
the outer regions of the circumnuclear disk are counter-rotating 
relative to stars in the inner disk and the ionized and molecu- 
lar gas. In this (putative) kinematically decoupled component 
(KDC), the stars in the outer parts are redshifted toward the NE, 
while the ionized+molecular gas and the stars in the inner parts 
are redshifted toward the SW. The separation between the two 
sets of kinematics is at roughly 7-8" radius, so the ring could 
be the separation between the tw o regimes. The most likely ex - 
planation for this, as proposed by lFalcon-Barroso et al.l ( 2006b . 
is that in NGC 5953 we are witnessing the formation of a de- 
coupled component as a result of an ongoing interaction; the 
gas in the inner parts could have been accreted from the com- 
panion NGC 5954. The KDC could be in its early stages, as 
judged by our stellar population age estimate of xL 10-15 Myr. 
This age would be consistent wit h the onset of the interaction 
of 40 Myr ago, as suggested by IJenkinsI (11984); this age from 
the time of closest approach is almost independent of the de- 
tails of the model. This kind o f counter-rotatio n has been also 
studied with simulations (e.g.. Idi Matteo et al.l 12008). and the 



Alternatively, the observed counter-rotation could be only 
apparent and not real, due to a warped disk in the center. The 
inner disk is not highly inclined, and a small tilt of the plane 
starting at a radius close to that of the ring, is able, because 
of projection effects, to reverse the apparent sense of rotation, 
without the presence of a true counter-rotation. A warped in- 
ner disk is expected from the close tidal interaction between 
NGC 5953 with NGC 5954, and the obvious tidal tails in op- 
tical and Hi. The warped disk hypothesis is supported by the 
fact that the ionised gas also, and not only the stars, shows 
some hint of counter-rotation. The existence of two gas flows 
counter-rotating in the same disk, is a very unlikely situa- 
tion, that could be at best transient. The gas counter-rotation 
just outside the ring i s visib le in the Sauron velocity maps of 
Falcon-B a rroso et al. (2006). Indeed, the velocity maps in HyS, 



and even more the [Oin] lines reveal on the NE part of the major 
axis that the velocity from blueshifted turns to redshifted again, 
by an amount of -40 kms . The apparent counter-rotation 
might look more striking on the stellar velocity maps, only be- 
cause the stellar component is less rotation-dominated in the 
inner disk, because of the higher velocity dispersion. So the 
velocities on the major axis turn from -40 km s" 1 to 40 km s _1 , 
while the gas is turning from -100 kms -1 to 40 kms -1 . But it 
is clear that in the NE part of the major axis, the stars and gas 
are co-rotating at the same velocities, and apparently counter- 
rotate relative to the inner disk. Since the gas cannot be counter- 
rotating with itself for a long time, the interpretation of a warp 
in the outer disk is probably more realistic. However, if the gas 
is actually counter-rotating at a radius of 1 .6 kpc, then it could 
arise from external accretion at that radius. This gas is config- 
ured in a full gaseous disk, of a bout 8 kpc radius, as mapped 
in Hi by IChengalur et al. ( 1994 ). We then can only witness a 
transient phase of settling of the accreted gas, which will soon 
align its rotation with the main disk, rotating in the same sense 
as the CO gas and stars in the center. 

5.4. Star formation in NGC 5953 

Most of the star formation in NGC 5953 occurs in the circum- 
nuclear region sampled by our 12 CO observations. The star- 
formation rate (SFR) given by Hff within a FOV of 33"x41" 
is O.6M yr _1 dFalcon-Barroso et al.l 12006). This corresponds 
to ~0.024M o yr~' kpc -2 . We measure a molecular mass of 
1 .7 x 10 9 M in a roughly similar region (see Sect. |4.2| >, which 
would give a molecular surface density of ~60M o pc~ 2 . In 
a Kennicutt-Schmidt (KS) law diagram without the neutral 
atomic gas component, this would place NGC 5953 in the 
transition regio n between norm al spiral disks and circumnu- 
clear starbursts dKennicuttll998l) . In fact, Hi observations show 
significant streaming mot ions over the inner 17 kpc inconsis - 
tent with global rotation dlono et al.l I2005E lHaan et al.ll2008l) . 
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Fig. 20. The strength Q (top) and phase d> (bottom) of the m 
= 1, 2, 3, 4 Fourier components of the potential, derived from 
the NIR image. The full lines correspond to m = 2 and the total 
strength, the dashed line to m = 1, dot-dash to m = 3, and dots 
to m = 4. 




offset along major axis (arcsec) 



Fig. 21. The 12 CO(1-0) contours are overlaid onto the gravi- 
tational torque map (t(x,y) x £(x,y), as defined in text) in the 
center of NGC 5953. The torque map (color scale) is plotted 
on a symmetric palette (wedge) although the maximum posi- 
tive torque is 50% higher in absolute value then the negative 
torque. The map is deprojected, and rotated so that the major 
axis of the galaxy is oriented parallel to the abscissa Ox. 



6. Gravitational torques on the molecular gas 

In this section we study whether gravitational torques, derived 
from the stellar potential in the inner region of NGC 5953, can 
account for the gas kinematics derived from CO and examine 
the efficiency of gravitational torques ex erted on the gas. A s 



descri bed in previous NUGA papers (e.g.. lGarcfa-Burillo et al 



2005), we assume that NIR images give the best approximation 



for the total stellar mass distribution, being less affected than 
optical images by dust extinction or stellar population bias. 



6.1. Evaluation of the gravitational potential 



Hence, the Hi probably does not take part in the star forma- 
tion activity in this region. It is probable that the enhanced star 
formation in NGC 5953 is caused by dynamical perturbations 
induced by the interaction. 

In NGC 5953, the high H2 surface density and the low 
SFR imply that the distribution of molecular gas (traced by 
CO) does not correlate locally with the SFR; the large reser- 
voir of molecular gas is not converted into stars proportion- 
ally to the KS law on small spatial scales. At present, the 
non-local correlation between SFR a nd gas is highly debated 
(Big iel et al.ll2008t iLerov et alj|2008l) . Other molecular transi- 



tions could be likely better indicators of SFR than the total H2 
content traced by CO. Tracers of dense molecular gas, such as 
HCN(l-O) and HCO + (1- 0) lines, are suspected to b etter cor- 
relate with the SFR (e.g. , |Gao & Solomonll2004alb1: IWu et al 
l2005HGracia-Carpio et alJ2008h . 



We computed the torques using Spitzer/IRAC and ARNICA 
images. They both yield comparable results, but we use in 
the following only the NOT/ARNICA images which give 
more spatial resolution. We perform the subtraction of fore- 
ground stars, deprojectio n, and resampling, as de scribed in 
other NUGA papers (e.g.. iGarcia-Burillo et alll2005l) . 

We repeat here some definitions and assumptions used to 
evaluate the gravitational torques. The NIR image is completed 
in the vertical dimension by assuming an isothermal plane 
model with a constant scale height, equal to ~ 1/1 2th of the 
radial scale-length of the image. The potential is derived by 
a Fourier transform method and assuming a constant mass-to- 
light (M/L) ratio whose value is chosen to better reproduce the 
observed 12 CO rotation curve. Beyond a radius of 20"(or 5.4 
kpc in diameter), the mass density is set to 0, thus suppress- 
ing any spurious m = 4 terms. This assumption is sufficient to 
compute the potential over the PdBI 12 CO(1-0) primary beam. 
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Fig. 23. The torque, or more precisely the fraction of the angular momentum transferred from/to the gas in one rotation -dhj L- is 
plotted for 12 CO(1-0) (left) and 12 CO(2-l) (right). The (red) dashed area corresponds to the resolution limit of our observations. 
In the left panel, the CO(1-0) map resolution is the limiting factor, while, on the right, the CO(2-l) and the similar NIR image 
resolution are the constraining factors. 
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Fig. 22. Same as Fig.|2T]for the 12 CO(2-l) emission taken as 
tracer of gas surface density. 



To characterize the non-axisymmetry of the potential, we 
expand its 0) in Fouri er components (m-modes), following 
Combes & Sanders! dl98ll) : 



*(/?, 0) = <S> (R) + V d> m (R) cos[m# 



m=l 



n(R)] 



where O m (/?) and <p,„{R) are the amplitude and phase of the m- 
mode, respectively. 

The strength of each m-Fourier component, Q m (R), is de- 
fined by the ratio between tangential and radial forces, Q m (R) = 
m<$> m /R\Fo(R)\. The strength of the total non-axisymmetric per- 
turbation is defined by: 

pmax(R) 

Qt(R) 



where Fj ax {R) and Fq{R) represent the maximum amplitude of 
the tangential force and the mean axisymmetric radial force, 
respectively. Figures l20l show the strengths and phases vs. ra- 
dius for the first m components. There is only a weak m = 2 
perturbation; the strongest one is m = I. 

6.2. Evaluation of gravity torques 

From the derivatives of <&(R,8) on each pixel we obtain the 
forces per unit mass (F x and F y ) and therefore the torques per 
unit mass t(x,y) can be computed by: 

t(x,y) = xF y -yF x . 

The torque map is oriented according to the sense of rotation in 
the plane of the galaxy, and the combination of the torque map 
and the gas density 2 map allows us to derive the net effect on 
the gas at each radius. The gravitational torque map weighted 
by the gas surface density t(x,y) x T,(x,y), normalized to its 
maximum value, is shown in Figures |2T] and [22] for 12 CO(l- 
0) and 12 CO(2-l), respectively. The observed gas distribution 
is representative of the time spent by a molecular cloud on a 
typical orbit at this location. 

To estimate the efficiency of the radial gas flow induced by 
the torques, we first compute the torque per unit mass averaged 
over azimuth, using "E(x, y) as the actual weighting function: 



t(R) 



J g nx,y)x(xF y -yF x ) 



F (R) 



where t(R) is, for definition, the time derivative of the spe- 
cific angular momentum L of the gas averaged azimuthally, 
t(R)-dhl dt \g. Now to have dimensionless quantities, we nor- 
malize this variation of angular momentum per unit time to the 
angular momentum at this radius and to the rotation period. We 
then estimate the efficiency of the gas flow as the average frac- 
tion of the gas specific angular momentum transferred in one 
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rotation (T rot ) by the stellar potential, as a function of radius: 
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where Lg is assumed to be well represented by its axisymmet- 
ric estimate, Lg = R x v ro/ . Figures [23] show AL/L curves for 
NGC 5953 derived from the 12 CO(1-0) (left) and the 12 CO(2- 
1) (right) data. The (red) dashed area corresponds to the res- 
olution limit of our observations. In the left panel, the CO(l- 

0) map resolution (~120pc) is the limiting factor, while, on 
the right, the CO(2-l) map resolution (~60pc) and the simi- 
lar NIR image resolution (~70pc) are the constraining factors. 
This figure shows that the torques are weak and predominantly 
positive between -100-400 pc in both 12 CO(1-0) and 12 CO(2- 

1) . Nothing can be said for the negative torques inside a radius 
of ~120pc for 12 CO(1-0) and ~70 pc, below our resolution. 
In 12 CO(2-l) the torques are noisier than in 12 CO(1-0), maybe 
because of a more clumpy emission. 

In summary, torques are predominantly positive in the re- 
gion ~1 00^4-00 pc for both lines, and the gas is not apparently 
fueling the central ~100pc region down to the effective spatial 
resolution of our observations. In addition, the absolute value 
of the torques is very small, less than 5% of the gas angular 
momentum is exchanged in each rotation. This is due to the 
very regular and almost axisymmetric total mass and gas dis- 
tributions in the center of the galaxy. 

7. Summary and conclusions 

The molecular gas in the Seyfert 2/LINER galaxy NGC 5953 
has been mapped with high resolution (2'.'1 x 1'.'4 for the 
12 CO(1-0) line and l'.'l x 077 for the 12 CO(2-l) line) in- 
side a radius of ~20" (~2.7kpc). The 12 CO emission is dis- 
tributed over a disk of diameter of ~16" with several peaks 
randomly distributed. The strongest intensity peak is not lo- 
cated in the center of the galaxy but shifted ~2-3" toward the 
west/southwest. The kinematics of the molecular gas show a 
general regularity with some local wiggles especially to the 
west of the nucleus maybe associated with the intensity peak 
toward the SW. 

By comparing the molecular gas distribution with obser- 
vations at other wavelengths, we studied correlations between 
different tracers of the ISM. Optical and NIR morphology has 
been studied by analyzing different images. The F606W/HST 
broad-band optical image of NGC 5953 shows a flocculent spi- 
ral structure and an "S" like-bar of ~250pc and i 60 pc in 
size, respectively. These structures are absent in the 1.6jum 
(HST/NICMOS/F160W) image and in the IRAC 8//m "dust- 
only" one; instead in the optical and at 1.6/zm we identi- 
fied at ~3" w est of center the foreground star studied by 
1990). The stellar structure has been investi- 
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gated by performing a two-dimensional B/D decomposition on 
the ground-based //-band image and IRAC 3.6 /im image. Both 
sets of residuals show a circumnuclear stellar ring with a radius 
of ~5-7", approximately coincident in size with the 12 CO disk. 
The presence of the ring in both sets of residuals suggests that 



it is not the typical "doughnut" associated with the diffraction 
limit and IRAC PRF incompatibilities, but rather is a real fea- 
ture. The size of this stellar ring is the same of the circumnu- 
clear star-forming ring seen in Ho- with previous observations. 
In addition, since this ring is visible in the NIR as stars but not 
as hot dust emission (J-K), we interpret it as a RSG population 
at least lOMyrold. 

There is an apparent counter-rotation between gas and stars 
inside the ring, and stars outside. The ring could be the sep- 
aration between the two kinematically decoupled components 
(KDC). The formation of a KDC could be explained as a re- 
sult of the ongoing interaction, and could be in its early stages. 
Alternatively, this apparent counter-rotation could be due to a 
warp of the plane of the disk. The ring would then be the start 
of the warp. 

Using NIR images we found that gravity torques acting on 
the gas are predominantly positive in the region ~1 00^-00 pc 
in both 12 CO(1-0) and I2 CO(2-l), indicating that the gas is not 
fueling the central ~100pc region, down to the effective spatial 
resolution of our observations. In addition, the absolute values 
of these torques are very small, certainly due to the almost ax- 
isymmetric total mass and gas distributions in the center of the 
galaxy. The AGN in NGC 5953 is apparently not being actively 
fueled in the current epoch. 
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